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ABSTRACT 

We report a correlation based on a spectral simulation study of the prompt 
emission spectra of gamma-ray bursts (GRBs) detected by the Swift Burst Alert 
Telescope (BAT). The correlation is between the -Epeak energy, which is the peak 
energy in the uFi, spectrum, and the photon index (F) derived from a simple 
power-law model. The Epcak - E relation, assuming the typical smoothly broken 
power-law spectrum of GRBs, is log^peak = 3.258 -0.829 F (1.3 < F < 2.3). We 
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take into account not only a range of -Epeak energies and fluences, but also distri- 
butions for both the low-energy photon index and the high-energy photon index 
in the smoothly broken power-law model. The distribution of burst durations 
in the BAT GRB sample is also included in the simulation. Our correlation is 
consistent with the index observed by BAT and E^^ak measured by the BAT, and 
by other GRB instruments. Since about 85% of GRBs observed by the BAT are 
acceptably fit with the simple power-law model because of the relatively narrow 
energy range of the BAT, this relationship can be used to estimate i?peak when 
it is located within the BAT energy range. 

Subject headings: gamma rays: bursts 



1. Introduction 



One of the fundamental characteristics of the prompt emission of gamma-ray bursts 
(GRB) is -Epeak; whlch is the peak energy in the vY^ spectrum. According to BeppoSAX and 
HETE-2 observations, . Epeak for GRBs is widely spread from a fe w keV to the MeV range 
as a single distribution (IKippen et al.ll2002l : ISakamoto et al.ll2005l ). This broad single .Epeak 
distribution strengthens the argument that these bursts arise from the same origin. Based 
on this observational evidence, there are several wo rks which try to understand a unified pic - 



ture of GRBs. For insta nce, the off-axis jet model (lYamazaki et al. 



the structured jet model ( iRossi et al. 



2002 



the variable jet opening angle model ( iLamb. Donaghy &: Grazianill2005l ) are the popular uni 



2004 



Toma et all 120051) 



Zhang fc Meszaro^ I2OO2I: Izhang et al.lbooi ). and 



fied jet models. On the other hand, there are theoretical models to explain the broad 



-E'peak distribution in the frame work of the internal shock model (IMeszaros et al 



Moch 


covitch et al.| 


2003 


I999I; 


Huang et al. 


2002 



2003I : iBarraud. et ahlbood ) and the external shock model (jPermer et al.l 



2002; Dermer and Mitman 



2002: 



20031 ) 



There are several important empirical relationships proposed based on the -Epeak energy. 
One of the most cited relationships is the correlation between E^^ak in the GRB rest frame 
{E^^ A and the isotrop i c radiated en ergy (i?iso), the so called the E^^^^-Eisq (Amati) rela- 



tion (Amati et al 



( ISakamoto et al. 



2002 



2004 



II 1 ^" ^ ' pi;a,JV — - V / 

Amatill2003l 1. Since this relation is extended down to X-ray flashes 
20061 ). the dynamic range of this relation is ~3 orders of magnitude 



in £^;'eak and 



The second correlation is between the -Epeak 



5 order of magnitude in E\. 
energ y and the collirnation- corrected energ y {E-,), the so ca l led E^ ^^^-E^ (Ghirlanda) rela- 

), this relation has much 



anda et al.l Em 



tion (IGhirlanda et al.ll2004l ). According to iGhir 
tighter correlation than the E^^^^^-E^^o relation 
ilar relationship, but without using E^ which is heavily dependent on the calculation of 



Liang fc Zhang (120051 ) investigated a sim- 
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the jet opening angle. They found a good correlation between -Epgak' -^iso, and the achro- 
matic break time in the afterglow light curve {tjet)- The third relationship is between -Epg^k 
and the isotropic peak luminosity [Lf^^^), the so called the E^e^k' -^fso^ (Yonetoku) relation 
( Yonetoku et al.ll2004l ). The latest fourth relationship is between Lf^^^, -Epeak^ the time 
scale of the brightest 45 per cent of the ba ckground subtracted counts in the light curve of 
the prompt emission (IFirmani et al.ll2006l ). If these relationships are valid, they must be 
related to the fundamental physics of GRBs. Thus, E^^^-^ energy provides us fruitful knowl- 
edge about the characteristics of the prompt emission of GRBs. Furthermore, knowing the 
-^peak energy is crucial to calculating the bolometric fiuence which reflects the total radiated 
energy in the prompt emission. 



A fter the launch o f Swi ft ( iGehrels et al I |200J) in 2004, the Burst Alert Telescope 
(BAT; iBarthelmy et al.l (120051 )) has observed about 100 GRBs per year. In about half of the 



GRBs , the Ep^ak energies are very likely to be within the BAT energy range (jSakamoto et al. 



2008al ). However, due to the relatively narrow energy band of the BAT (15-150 keV in the 
background subtracted spectrum using the mask modulation), the BAT has a difficulty in 
determining -Epcak- purpose of this study is to find a way to estimate E°]^l^ when it lies 
within the BAT energy range. 

Here, we report a good correlation between the photon power-law index derived from a 
simple power-law model and i?poak based on the spectral simulation study. We use a sample 
of 31 long BAT GRBs that are well fitted with the power-law times exponential cutoff model, 
and also 26 GRBs observed by other GRB instruments concurrent with the BAT to confirm 
our correlation. Our correlation provides an estimate for -Epeak from the photon index in a 
simple power-law fit at the range from 1.3 to 2.3. We also calculated the la confidence level 
of the estimated i?peak of our correlation. 



2. BAT spectral simulation 



Because of the syst ematic difference i n the spectral parameters based on the assumption 
of the spectral model (IBand et al.l Il993l ). we decided to perform the simulations for two 
typical G RB spectral n iodels as input spectra: the smoothly broken power-law model (Band 
function; band et al.l Jl993h ) and a power-law times exponential cutoff modej^ (GPL) 
model. We fit the low-energy photon index, a, and high-energy photon index, (3, of 124 



MN/dE = KiE^i cxp[-E(2 + ri)/Epeak] if E < (Ei 

(ri-E2)-Bpeak/(2 + ri) 

MN/dE ^ E" cxp(-(2 + a) E/Ep^ak) 



r2)£;poak/(2 + El) and dN/dE 
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sampl es of the Band function fit ("BAND" in their notation) in Table 9 of iKaneko et al. 
( 120061 ) by the normal distribution. We obtained a of —0.87 with o of 0.33 and of —2.36 
with (J of 0.31. N ote that we are not excluding the case of /3 > —2 in our simulations because 
two reports (e.g., ISato et al.ll2005l : iKaneko et al.ll2006l ) show fits with (i > —I'm. both time- 
averaged and time-resolved burst spectra. However, the fraction of simulated spectra with 
/3 > —2 is only 13% of the total. Similar ly, for a CPL model, we fit the low-energy photon 
index, acPL, for the sample in Table 9 of IKaneko et al.l (120061 ) ("COMP" in their notation; 
67 samples) by the normal distribution. We found «cpl of —1.11 with o of 0.30 (see Figure 
[1]). These a, /?, and acpL distributions are used in our spectral simulation. 

In our simulations, -Epeak varies from 1.4 keV to 1210 keV in a logarithmic scale. The 
15-150 keV fluence varies from 5 x 10~* to 5 x 10~^ ergs cm~^ in a logarit hmic scale. The 



fluence range is determined based on the BAT observations (BATl catalog; ISakamoto et al. 



(l2008al )). The simulation used 20 values for fluence and 70 values for .Epcak- The exposure 
time of the spectrum is the best fit log-normal distribution of the BAT Tioo duratioiifl 
reported in the BATl cataloj^ (See the bottom panel of Figure [1]). The normalization of 
the input spectrum is calculated to be the input fluence value. The spectral simulations 
are performed 1000 times for each grid point. The background is included in the simulation 
using the spectrum created from the event data of the false BAT trigger 180931. Since the 
background is subtracted using the mask modulation, the exposure time of the background 
spectrum is set as the same as the duration of the foreground spectrum. Four incident angles, 
on-axis (0°), 15°, 30°, and 50° off-axis, are simulated independently. The simulated spectra 
are fitted from 14 keV to 150 keV with a simple power-law models! (PL), a CPL, and the 
Band function. Xspec 11.3.2 was used in both creating and fitting the simulated spectra. 



Figure [2] shows the numbers of the simulated spectra which have A^^ ( A^^ 



2 

PL 



xland Band function or A^^ = Xpl ~ Xcpl ^ CPL model) greater than co as a 

function of the -Epeak and the energy flux in the 15-150 keV band. This A^^ > 6 corresponds 
to >2.4 a confidence. The figures in the left and right row show the results based on 
the Band function and a CPL model, respectively. We note the distinct differences in the 
shapes of the confidence contours, especially at low .Epeak, between the Band function and a 
CPL model as an input spectrum. The results show that if a CPL model is indeed a true 



•^The duration includes from to 100% of the GRB fluence. 
"^The duration between tstart and tstop time of the fluence table. 
5dN/dE - E"^ 

^This is a current criterion used in the BAT team for reporting the spectral parameters based on a CPL 
fit in the BAT refined circular of the Gamma-ray Burst Coordinates Network. 
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spectral shape, BAT can measure Epcak at the lower boundary of its energy range (~ 15 keV) 
with a very high significance. On the other hand, a low .Epeak measurement would be very 
challenging if the Band function is the true spectral shape. Figure [3] explains the reason for 
these differences. The figure shows the calculated photon spectra in a GPL model and the 
Band function for E'pcak = 15 keV. In a GPL model, the spectrum can not be fit with a PL 
model because of the curved shape (exponential component) in the BAT observed energy 
band. Therefore, we would expect a significant improvement in with a GPL fit over a 
PL fit. However, in the Band function, due to the extra power-law component (high energy 
power-law component) in the formula, the spectrum at the BAT observed energy band would 
be just a simple power-law with a high energy photon index. This is the reason why we see 
a difference in the confidence contours based on the assumed spectral models. The results 
also show that the .Epeak measurement becomes difficult for BAT when i?pcak is below 30 keV 
or above 100 keV in the Band function shape. In the GPL shape, .Epeak can be determined 
even at ~ 15 keV. 

Figure H] shows the number of the BAT GRBs which can be acceptably fit by a PL 
model and by a GPL model as a function of the 15-150 keV fiuence. The data are from the 
BATl catalog. In the case of an incident angle less than 25 degrees, a GPL model becomes 
an acceptable fit for fiuence >10~^ ergs cm~^. However, a PL model still be acceptable fit if 
Epcak is located above or below the BAT energy range. On the other hand, the fiuence must 
typically be greater than 3 x 10~^ ergs cm~^ in the case of an incident angle greater than 50 
degrees. These threshold fluences required to measure Ep^ak in the BAT data correspond to 
the ~50% confidence contour (green) in our simulation results of Figure [2J 

Next, we made histograms of Epc^k for each photon index on a 0.1 grid from to 3.5 
using the range of fiuences corresponding to the 1-cr interval of the BAT observed fiuence 
distribution in the BATl catalog. The 1-cr fiuence interval corresponds to the range from 
3.4 X 10"^ ergs cm~^ to 5.4 x 10"^ ergs cm~^. This selection of the fiuence range allows us 
to reduce the systematic effect of the inclusion of unrealistically bright or dim simulations. 
Furthermore, since we are interested in estimating the .Epeak for the bursts which do not 
show a significant improvement in by a GPL fit over a PL fit, we also only selected the 
simulated spectra with Ax^ = Xpl~Xcpl < 6. Because the numbers of simulated spectra are 
different for each i?pcak grid due to these selections, we normalized the number of simulated 
spectra in each P-i^peak grid by the total number of spectra in each i?peak grid. Figure [5] 
shows the contour map of the photon index (P) and logiJpcak for the Band function (left) 
and a GPL (right) model. There is a correlation between P and log£^peak in the range from 
1.3 to 2.3 of P for the Band function. The correlation continues to P = 3.0 in the case of a 
GPL model for the same reason as we demonstrated in Figure [31 It might be interesting to 
note that a very steep photon index such as P ~ 3 is not possible to achieve if the source 
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spectrum is the Band function. In this case, the source spectrum might be much closer to 
a CPL shape. One important conclusion is that the correlation between logE'pcak and F 
exists independent of the incident angle of the burst. Therefore, this correlation, the -Epeak 

- r relation, can be used for all BAT long GRBs within the allowed F range, although with 
larger uncertainty for GRBs at large incidence angle. 

We extracted the peak i?peak value from each histogram of F and fit with a linear function 
using the range from 1.3 < F < 2.3 for both the Band function and a CPL model. Although 
the correlation exists until F = 3 in a CPL case, we use the same F range for the Band 
function and a CPL model to investigate the systematic difference based on the assumption 
of the source spectrum. The best fit i?peak-F relations are summarized in Table [1] (Band 
function) and Table [2] (CPL). To estimate the 1-a uncertainty of the relation, we found 16% 
and 84% points of -Epcak from each histogram of F and fitted with a cubic function from 
1.3 < F < 2.3. The best fit cubic functions of the lower and higher 1-a confidence level are 
also summarized in Table [1] and Table [2l Figure [6] shows the best fit functions of the i?peak - 
F relation and its 1-a confidence level with the data points used in the fittings. We note that 
the wide i?peak range in the simulations (in our case from 1.4 keV to 1210 keV) is essential 
to derive the 1-a confidence level of the relation. If the -Epeak range in the simulations is 
not wide enough such as from 10 keV to 500 keV, we noticed that the confidence level will 
be underestimated by a factor of 2 for the upper limit at F of 1.3 and by a factor of 5 for 
the lower limit at F of 2.0. Due to the smoothly curved shapes of the Band function and 
the CPL model, the -Epeak grids in the simulations have to be an order of magnitude wider 
than the energy range of the instrument, so that a curvature (or i?pcak) in the spectrum 
is completely outside the energy range of the instrument for the i?pcak around the energy 
limits of the instrument. However, we also notice that the best fit ii'peak - F relation itself 
is less sensitive to the energy limits on the simulations. Although the confidence level is 
different between the Band function and a CPL model, the best fit linear function shows 
little difference between these two spectral models. We also calculated the relation weighting 
the results at the incident angles of 0°, 15°, 30°, and 50° by the distribution of the incident 
angle of the BAT GRBs (Figure [7]). Hereafter, we call this relation as the weighted -Epeak 

- F relation. The contour plots of the weighted -Epcak - T relation, the plot of the best fit 
functions, and the formula of the best fit functions are shown and summarized in Figure [HI 
Figure [H and Table [1] and Table [21 respectively. 

In the application of our -Epeak - F relation, we strongly encourage the reader to use 
the result based on the Band function as a prior. The main reason for also performing 
the simulations of a CPL model as a prior is to see the systematic effect due to a prior 
assumption of the spectral model. Our results are clearly demonstrating the effect of the 
assumed spectral model. From the various measurements of the burst spectra by different 
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instruments, the true burst spectrum is very likely to be the Band function at least for long 
GRBs. Therefore, the -Epeak - T relation based on the Band function as a prior is the most 
suitable relation to apply for the BAT long GRBs. 



3. Comparison to other i?pcak Measurements 



To investigate the validity of our sirn ulation study, we used the spectral parameters on 
the BATl catalog (jSakamoto et al.ll2008al ). Table [3] shows the spectral parameters of 31 long 
GRBs (Tgo > 2 seconds) having Ax^ of greater than 6 in a GPL model over a PL model fit. 
Figure [TD] shows Ep^ak energy in a GPL model and the photon index, F, in a PL model for 
the BAT GRBs overlaid with the weighted .Epeak - T relation. We also plot F derived from 
the BAT data and Epcak reported by Konus-Wind or HETE-2 in the Gamma ray bursts 
Goordinates Network (GGN) listed in table HI As seen in the figure, the 1-a confidence level 
of the -Epeak - r relation based on the simulation study is consistent with the 90% confidence 
level of F observed by the BAT and -Epcak observed by the GRB instruments. However, we 
want to caution about using our i?peak - F relation for estimating i?peak- Our estimator is 
based on prior assumptions of the low-energy and/or the high-energy photon index measured 
by the BATSE. Therefore, -Epeak based on our i?pcak - T relation only provides a likelihood 
of the -Epcak value not the "measurement." 

The calculation of the bolometric flux or fiuence is another challenge when using the 
BAT data alone. However, since the low energy photon index a and the high energy photon 
index (3 of the Band fu nction are quite stable parameters even if -Epcak varies from a few 



keV to a few MeV (e.g., ISakamoto et al.ll2005l : iKaneko et al.ll2006l ). one could estimate the 



bolometric flux or fluence assuming the best fit a and (3 from the BATSE time-averaged 
spectral analysis of iKaneko et al.l (120061 ) . and applying the best fit -Epcak derived from our 
-Epeak - F relation. To get the normalization for the Band function spectrum, one would scale 
the Band spectrum so that the flux in the BAT energy range matches to the BAT measured 
flux. One can also estimate the error of the flux or fiuence by propagating the errors of a, 
(3, -Epcak, and the normalization, however this estimate will not be strictly correct because 
the parameters of the Band function are correlated. Finally, we caution against relying too 
heavily on this derived bolometric flux or fiuence, since the method described uses averaged 
a and (3 from a different burst population and an estimated, rather than measured -Epcak- 
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Discussion 



According to ISakamoto et al.l (120051 ). an equal number of X-ray flashes (XRF), X-ray- 
rich GRBs (XRR), and GRBs are reported in the HETE-2 GRB sample. Their classiflcation 
of GRBs based on the fluence ratio between the 2-30 keV and 30-400 keV bands is almost the 
equivalent of classifying GRBs by Ep^ak- The boundaries of .Epeak between an XRF and an 
XRR, and an XRR and a GRB are around 30 keV and 100 keV. When we use the weighted 
-E'pcak - r relation for the Band function to calculate the corresponding F for each Epeak, T 
is ~ 2.2 and ~ 1.5 for Epcak of 30 keV and 100 keV, respectively. Applying these F criteria 
to a sample of 206 BAT GRBs, excluding short GRBs (Tgo < 2 seconds) and GRBs with 
incomplete dataset, we found that the number of XRFs, XRRs, and GRBs are 20, 126, and 
60 respectively. The numbers of XRFs, XRRs, and GRBs in the HETE-2 sample are 16, 
19, and 10 respectively. Therefore, the ratio of the numbers of XRRs and GRBs is identical 
for both BAT and HETE-2 sample. The small numbers of XRFs in the BAT sample is due 
to the difficul ty in o b servin g very soft XRFs in the BAT (IBand I l2003l . |2006| ). However, as 
mentioned in iBand I (120061 ). it is very difficult to determine the actual detection threshold 
of the BAT due to its complexity in the triggering algorithm. Although nothing could be 
addressed about the actual number of XRFs, the number of GRBs in XRRs and GRBs seen 
in the BAT sample is consistent with the HETE-2 sample. The detailed study of the Swift 
XRFs and XRRs is presented elsewhere (ISakamoto et al.ll2008bl ). 



Butler et al.l (120071 ) calculated -Epeak by their Bayesian approach for 218 Swift GRBs 



using only the BAT data. Based on their calculated Epp^ k and the bolometr ic fluence, they 
claimed that all of th e empirical relations, E^f^y - E\^n (Amati et al.ll2002l ). E^^^^^ - Lf^f^ 
( Yonetoku et al.ll2004l ). and E^^^^T^r, - Lf^^^ ( Firmani et al.ll2006l ). proposed in the pre- Swift 
observations are not valid for the Swift BAT sample. We investigated the validity of their 
Epeak by checking i?peak obtained by using our i?peak - T relation. \ ye created th e BAT spectra 
for their GRB samples by the time interval reported on Table 1 of lButler et al.l (120071 ). Then, 
we flt the spectrum by a PL model to extract the best flt F. By only selecting their F within 
the allowed F range for applying our i?poak - T relation (1.3 < F < 2.3) and also excluding 
the short GRBs (156 samples in total), we calculated -Epcak applying ou r weighted Epp^k - F 
relation for the Band function. Figure [11] shows the -Epeak reported on iButler et al.l (120071 ) 
versus -Epeak derived from our weighted -Epeak - T relation for the Band functio n. Although 



the er ror bars are large in both estimators, the flgure shows that -Epeak of the iButler et al. 



(120071 ) sample has a systematica l ly hig her -Epeak compared to that from our -Epeak - T relation. 
About 20% of the IButler et al.l (120071 ) sample selected based on the range of F from 1.3 to 
2.3 exceeds -Epeak ~ 150 keV which is the limit of the estimated -Epeak using our -Ep e ak - F 
relation for F = 1.3. Furthermore, we are already excluding 20% of the iButler et al.l (120071 ) 



sample because those bursts fall outside limit range of F from 1.3 to 2.3 in our relation. This 
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limit is determined because -Epeak is very likely located outside of the BAT energy range, and 
therefore, the BAT data alone can not constrain about .Epsak (the BA T data only can provide 
the limit in i?peak)- In total, about 35% of the iButler et all (120071 ) samples are obviously 
inconsistent with the i?peak es timated based on o ur i?pcak - T relation. However, -Bpeak is 
constrained in the most of the iButler et al.l (120071) sample . These results provide a caution 
for the method for estimating i?peak in iButler et al.l (120071 ). 



Butler et al.l (120071 ) justify their Swift-only Ep^ak estimates in part by comparing to 



Konus- Wind measurements of i?peak for the same bursts. They also using the .Epeak distri- 
bution as a prior. However, the assumpt i on of -Epeak measured by Konus- Wind should be 
identical to that of BAT in iButler et al.l (120071 ) might not be valid. Because BAT has a 
significantly larger effective area and also relatively softer energy band than Konus- Wind, 
the time interval for creating the time-averaged spectrum based on the BAT data could 
be systematically longer than that of Konus- Wind (Sakamoto et al. in preparation). This 
longer time interval for the time-averaged spectrum in BAT might lead to a systematically 
lower -Epeak which might contradict with the .Epcak based on the Konus- Wind data alone. 
For instance, GRB 060117, which has individual measurements of i?peak from the BAT and 
the Konus- Wind data (see Table 2 and Table 3), shows a smaller i?poak in the BAT data. 
The d uration reported based on the Konus- Wind data is ~ 20 seconds (iGolenetskii et al. 



2006bl ). On the other hand, th e duration used to accumulate the BAT spectrum is ~ 30 



seconds (ISakamoto et al.ll2008al ). We confirmed based on our cross-calibration work that 
there is no systematic difference in Ep^^k of this burst between BAT and Konus- Wind if we 
select exactly the same time interval for accumulating the spectrum (Sakamoto et al. in 
preparation). Therefore, we believe that a prior assumption of E'pcak based on a particular 
GRB instrument might introduce an another level of a systematic error in the analysis. Most 
importantly, we believe that testing these empirical relations, which require the broad-band 
spectral properties of the prompt GRB emission, by using only the BAT narrow-band data 
could lead to a wrong conclusion. Current on-going activity for analyzing the spectral data 
of simultaneously observed BAT GRBs by other GRB missions such as Konus- Wind and 
Suzaku/WAM (Sakamoto et al. in preparation; Krimm et al. in preparation) is indeed a 
necessary step to answer for the validation of these empirical relations. We might want to 
emphasize that the 1-a confidence level of our i?pcak — T relation based on the Band function 
includes most of i?pcak reported by other instruments (see Figure [TOj) . Therefore, the confi- 
dence level which we are quoting in our estimator is large enough to include the systematic 
problem in -Epeak among the different instruments. 

We report the correlation between i?peak and the photon index, F, of the BAT prompt 
emission spectrum based on our simulation study. Using this relation, it is possible to 
estimate i?peak from F in the range from 1.3 to 2.3. We also performed the spectral simulations 
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for assuming various incident angles (0°, 15°, 30° and 50°) and different spectral models 
(Band function and CPL). However, none of these systematic effects changes the relation. 
In the application, the i?pcak - T relation based on the Band function as a prior is the 
appropriate formula to use. The -Epeak - T relation could be informative for classifying the 
BAT GRBs from the photon index alone as derived from a simple power-law model which 
is the best fit for about 80 % of the whole population of the BAT GRBs. 



We would like to thank the anonymous referee for comments and suggestions that ma- 
terially improved the paper. 
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Table 1. i?pcak - T relation based on the Band function 
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Table 2. i?pcak - T relation based on a CPL model 
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Table 3: The BAT time-averaged spectral parameters fitted with a simple power-law (PL) 
model and a power-law times exponential cutoff (CPL) model. See the BATl catalog paper 
for the details about the BAT analysis (Sakamoto et al. 2008a). The degree of freedom in a 
PL fit and a CPL fit is all 57 and 56 respectively. 
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"The confidence interval is not calculated because of > 2. 
Short GRBs, GRB 050820A and GRB 050925 are excluded. 
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Table 4: The spectral parameters of simultaneously observed by Konus- Wind or HETE-2. 
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Fig. 1. — Input parameters in the spectral simulations. The distribution of the low energy 
photon index a, the high energy photon index /3 in the Band function, the low energy photon 
index a in a CPL model from the BATSE GRB sample, and the BAT Tioo duration from 
top to bottom, respectively. The dotted line represents the best fit in a gaussian. 
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Fig. 2. — The contour maps showing, as a function of the simulated -Epeak and the energy 
flux in the 15-150 keV band, the number of simulated spectra which have Ax^ > 6 (left row 
for the Band function: Ax^ = Xpl " XBandi right row for a CPL model: Ax^ = Xpl ~ 
The incident angles of the simulations are 0°, 15°, 30° and 50° from top to bottom. 
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Fig. 3. — The schematic drawing of the photon spectra of the Band function (black) and a 
CPL model (red) with .Epeak of 15 keV. The low energy photon index is —1 for both models. 
The high energy photon index of the Band function is —2.3. A PL model with a photon 
index of —2.3 is also overlaid in the plot (blue). The vertical dotted lines are the BAT 
observed energy band of 15-150 keV. 
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Fig. 4. — The number of GRBs acceptably fit by a PL model (top panel) and by a CPL 
model (bottom panel) as a function of the fiuence in the 15-150 keV band. The samples of 
the incident angles of bursts less than 25 degrees (top) and larger than 50 degrees (bottom). 
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Fig. 5. — Contour maps showing the number of simulated spectra as a function of photon 
index and input i?peak- The left and right rows are the Band function and a CPL model, 
respectively. The incident angles of the simulations are 0°, 15°, 30°, and 50° from top to 
bottom. 
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Fig. 6. — The best fit Epcak - T relations (solid line) and the lower and higher 1-a confidence 
level of the relations (dashed lines) for the Band function (top) and a CPL model (bottom) 
with the data points (circles: -Epcak at the peak of the histogram of F, squares: i?poak value 
of 16% crossing point of the histogram of F, and triangles: i?peak value of 84% crossing point 
of the histogram of F). The black, red, green and blue show the cases of incident angles 0°, 
15°, 30° and 50°, respectively. 
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Fig. 7. — The incident angle (6) distribution of the BAT GRBs. The dotted hne is the best 
fit gaussian model. 
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Fig. 8. — Contour maps showing the number of simulated spectra as a function of photon 
index and input .Epeak after weighting the simulation results of 0°, 15°, 30°, and 50° incident 
angles by the incident angle distribution of the BAT GRBs shown in Figure [7] (top: the Band 
function and bottom: a CPL model). 
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Fig. 9. — The best fit weighted .Epeak - T relations by the incident angles (solid line) and the 
lower and higher 1-cr confidence level of the relations (dashed lines) for the Band function 
(top) and a CPL model (bottom) with the data points (circles: i?poak at the peak of the 
histogram of F, squares: E'peak value of 16% crossing point of the histogram of F, and 
triangles: -Epeak value of 84% crossing point of the histogram of F). 
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Fig. 10. — The distribution of -Epcak and photon index, F, in a PL fit in the BAT GRB sample 
(black circles). The GRBs which were simultaneously observed by HETE-2 and Konus-Wind 
are overlaid (blue triangles). The weighted -Epcak - T relation for the Band function (top) 
and a GPL model (bottom) with 1-a confidence level is overlaid on the data. Inset: The 
extended figures of F from 1.3 to 2.3 where the -Epeak - F relation is valid. 
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Fig. 11. — The relationship between .Epeak reported by lButler et al.l (120071 ) and Epcak derived 
from the weighted Epcak - T relation for the Band function. The sample only contains 
lor ig bursts wh i ch hav e a PL photon index F from 1.3 to 2.3. Both i?peak distributions 
of iButler et al.l (120071 ) and the -Epeak - T relation for the same sample are shown in the 
histograms. 



